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Abstract. Compressed sensing investigates the recovery of sparse signals from linear measurements. But 
often, in a wide range of applications, one is given only the absolute values (squared) of the linear mea¬ 
surements. Recovering such signals (not necessarily sparse) is known as the phase retrieval problem. We 
consider this problem in the case when the measurements are time-frequency shifts of a suitably chosen 
generator, i.e. coming from a Gabor frame. We prove an easily checkable injectivity condition for recovery 
of any signal from all time-frequency shifts, and for recovery of sparse signals, when only some of those 
measurements are given. 
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1. Introduction 

Phase retrieval, a common problem in a variety of applications including X-Ray crystallography, optical 
imaging and electron microscopy, is the task of recovering a signal from the squares of the absolute values of 
its linear measurements. The best one can hope for in this case is to recover the signal up to a unimodular 
constant, because x and cx, where |c| = 1 will always give the same measurements. To fix notation, let 
P = ^ bs a set of measurement vectors, where K is K. or C. Further, let T = {c G K : |c| = 1}. 

The measurement process is then given by the map 

MFix)=[\{x,h)\^ |(a;,/2)|2 ... \{xjm)\^f. 

The task is to recover a: up to a global phase, given j\4f(x). We say that F allows phase retrieval, if the 
map MF is injective. 

There are three main directions of questions that one is interested in when looking into this problem: 

• Injectivity: Which properties of the measurement vectors can give us necessary and/or sufficient 
conditions on the injectivity of the map Mf? 

• Minimal number of measurements: How many measurements are needed for a set F to allow phase 
retrieval? 

• Algorithms: How can one practically find x, given the intensity measurements A4f(x)? 
Comprehensive answers to these questions and open problems can be found in gUSUlls]. Motivated by 
different applications, these general questions can also be asked for only a particular type of measurements, 
and/or signals. In this paper we focus on Gabor, or differently said, short-time Fourier measurements, which 
are time-frequency shifts of a suitably chosen generator. This type of measurements is of particular interest 
for many applications in speech and audio processing EH, ptychographical GDI m etc. On the signal side, 
we are looking at the sparsity constraint - a natural assumption that the signal we want to recover has 
only few non-zero entries, or is a linear combination of few vectors (has sparse representation). The sparsity 
constraint is a novel paradigm for signal- and image processing, and utilized, in particular, for compressed 
sensing methodologies m- Sparse phase retrieval is studied in EH ES] ■ 

A combination of phase retrieval from Gabor measurements for sparse signals was firstly considered in 
m- The theoretical results there are about the recovery of non-vanishing signals from a full set of Gabor 
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measurements, and some intuition about the difficulty of recovery of sparse signals is given. Numerical results 
show that recovering sparse signals can be effectively conducted with modification of the GESPAR algorithm 
[21] . using less than N'^ measurements. 

In the very recent work m both theoretical and numerical investigations show, that 0(A^log^(7V)) 
measurements are enough for recovering general signals from block circulant Fourier based measurements, 
and if the signal is A:-sparse, only 0(A: log^(A^)) measurements. The structure of the measurements is similar 
to the one of Gabor systems, but at this moment it is not clear how their results transfer to the Gabor 
setting that we consider here. 

In this paper, we investigate the problem of phase retrieval from Gabor measurements, for both full 
and sparse signals. Our main concern is the question of injectivity. Using the characterization of phase 
retrievability via the properties of the kernel of the PhaseLift operator [3, we provide a condition on the 
generator, sufficient for the corresponding Gabor system to allow phase retrieval. We show how this condition 
can be eased, if the signal that needs to be recovered is non-vanishing. We further provide two examplary 
classes of generators, complex random signals and characteristic functions of difference sets, which satisfy the 
above mentioned condition. The common Gabor generators, which are short windows or Alltop sequences, 
on the other hand - as we will show - are not suitable for phase retrieval of general signals (they fail to 
recover sparse signals), this problem was also considered in |17] . 

Further, we extend the injectivity condition from to the sparse setting, and provide a similar, but 
more involved condition on the generator which can guarantee us phase retrievability of sparse signals, 
additionally with less than measurements. We generalize this result also to signals which are sparse in 
a Fourier domain. When N is prime, we construct generators such that the Gabor system can do fc-sparse 
phase retrieval from 0{k^) measurements. As we will see, the result can also be interpreted as an injectivity 
condition for recovery of structured fc^-sparse vectors from 0{k^) linear measurements. 

Both the injectivity theorems naturally provide a simple algorithm for recovery of signals up to a global 
phase. When all N'^ measurements are given, the recovery of any signal is possible by using solely the 
fast Fourier transform, making the algorithm extremely fast. If some of the measurements are lost, we can 
employ £i minimization to get the signal back. We provide several numerical experiments to test this idea 
for various settings. Although the number of measurements for recovery in our work is still of relatively high 
order, the ideas we use are novel, and might be of interest for the community for better understanding of 
this problem, and its future development. 

The remainder of this paper is organized as follows: Section[2]is dedicated to the phase retrievability ques¬ 
tion for general signals, from all Gabor measurements. In Section [31 we focus on the sparse setting, and 
show that /c-sparse phase retrieval is possible with order of k^ Gabor measurements. A detailed description 
of the algorithm that we propose, and its empirical evaluation is presented in Sectional 


2. An injectivity condition for arbitrary signals 

2.1. Notation and basic objects. We will be working in the signal space C^, as a space of complex 
valued, N periodic functions with integer argument, x = x{j), j S Z, which therefore always has to be 
assumed modulo N. We will use the customary domain [0,..., — 1] of j, but we will often write j G Zjv 

for convenience. The scalar product between two signals x and y is defined as 


N-l 

{x,y) = X! 

j=o 

and the Hilbert Schmidt scalar product between two N x N matrices A and B as 

=ti{A*B) = ^ (Aei.Bei). 
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The A^-th root of unity will be denoted by w = e^. We define the (discrete) Fourier transform x and the 
inverse Fourier transform x oi x G as follows: 

N-l ^ N-1 

H x{j) = 

n—0 n—0 

A family of vectors in is called a finite frame for [7], if there exist constants 0 < A < i? < oo 

such that 

M 

A||a:||^ < ^ \{x, < -B||a:||^ for all x G C^. 

i=l 

If A = B is possible, then is called an A-tight frame. 

For p G In, we define the translation operator Tp : i-A through 

{Tpx){n) = x{n-p) 

Further, we define for £ G Ijn, the modulation operator Mi : i—>■ C'^ through 

{Mix){n) = u}^'^x{n). 

A Gabor frame 0 is the collection of all translations and modulations of a single vector g G C^, 

iMiTpg)^-Jo- 

For a pair A = (jp,£) sometimes we will use the short-hand notation !!> := MiTp, and g\ := IlAg. The 
matrices of those operators are unitary, and the collection of them forms a basis of [23], i.e. 

{n^,Il^)p^,=NSp,N ( 2 . 1 ) 

We will need the following well-known commutation relations between translations and modulations. 
Lemma 2.1. (TS) Let X = {p,l),pL = {q,j) G Z^. Then, we have 

MeTp=oj^PTpMi, 

2.2. Injectivity for full Gabor measurements. As mentioned before, we want to pose the question 
under what conditions a signal x from some class C C can be recovered from a set of its Gabor intensity 
measurements (|(a:, (/A)|^)AeA, A C Z^. Since these measurements are invariant under multiplication with 
c G T = {c G C, |c| = 1}, the best we can hope for is to recover x up to a global phase. If we denote by C-^/T 
the set of equivalence classes under the equivalence relation x^y<t^3cGT:x = cy, we can formally pose 
the problem as follows: Under what conditions on g is the map 

7Wg:C/T^R^', (|(x,gA)|") AgA 

injective? 

Defiuitiou 2.1. We say that the Gabor system G = (5'A)AgA associated to a generator g G is allowing 
phase retrieval for C (or has the phase retrieval property), if the map AIg is injective. 

We start by considering the problem of recovering arbitrary signals from all measurements, i.e. C = 
and A = Z^. In order to investigate which Gabor frames are allowing phase retrieval for this class, we 
will use a well known characterization of the phase retrieval property in the complex case, given via the 
properties of the kernel of the PhaseLift operator, also called super analysis operator in |^. For a set of 
measurement vectors {fi)^i in this operator is defined as 

: C^xAT ^ C"*, H^{{H,fif:)j,s)T=i^ (2-2) 

Note that when H is in the form xx*, {H, fifi)NS ~ ^fi) — /i)P- note that the authors of [5] 

chose the set of Hermitian matrices as the domain of A. We define A in this way in order to avoid some 


^We in fact always obtain a frame in this way if 3 7 ^ 0. The frame is even N || 3 ||^-tight. m 
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technicalities. However, the space of Hermitian matrices is a very natural domain in the context of phase 
retrieval, as the next theorem suggests. 

Theorem 2.1. [S] A set of measurement vectors allows phase retrieval if and only if the kernel of 

the associated map A does not contain any Hermitian matrices of rank 1 or 2. 

With this theorem, we can prove that the full set of Gabor intensity measurements allows phase 
retrieval, as long as a simple condition is satisfied. 

Theorem 2.2. Let g G be a generator for which 

(2.3) 

for every X G Then the corresponding Gabor frame G = {gx)x^TP allows phase retrieval. 

Proof. Theorem 12.11 suggests that we should investigate {g\,Hg\) for H G Equality (12.11) implies 

that 

^ ~ ^ n^. 

If // = (p, ^), we have 

where Hp denotes the (discrete) Fourier transform of the vector Hp, defined by 'Hp{i) = Hi^i-p. Note that 
Tip is in some sense the p-th ’band’ of the matrix H. It hence holds 

N {H, gxgDns = N {gx, Hg^) = ^ (gx^vm^gx) = E ^ ^^9^) ■ 

p,£ p,£ 

If we write A = {q,j)) we know by Lemma [2.11 that Il(^p t^^g\ = n(p £)n(q j)n(p Using 

this, and the fact that n>, is unitary, we arrive at 

N {gx, Hgy) = ^ {g, gp.i) ■ (2.4) 

p,£ 

Now, assume that this vanishes for all A = (<?, j) G Fixing j, we see that the above expression is just 
the value of the Fourier transform of the vector G with pth entry 

U«(p)=E^''^pW (9,9p.i) 

i 

evaluated at j. Since dH equals zero for all j, the vector U® vanishes for every q. Further, we observe that 
V‘^{p) is N times the value at q of the inverse Fourier transform of the vector G , where 

wP{l)=np{e){g,gp,i). (2.5) 

This expression must therefore be equal to zero for all p and £. With the assumption on the generator, we 
conclude that all the vectors Hp must vanish, and therefore also H. H can hence not have rank 1 or 2, and 
the proof is finished. □ 

Carefully going through the argument of the last proof, we see that it shows that the only matrix in the 
kernel of A is the zero matrix. Therefore, the proof actually shows that, under the assumption (12.3L A is 
an injective map. We use this idea to prove the following theorem. 

Theorem 2.3. Let g G be such that {g,gx) 7 ^ 0 for all A G Then, the rank-1 operators 

{9xgx)x^iP form a frame for ^^ (equipped with the Hilbert-Schmidt norm) and hence a basis. The frame 
bounds are given by 

A = N ■ min |(5,5A)|^ B = N ■ ma.^\{g,gx)\^ . 


PHASE RETRIEVAL FROM GABOR MEASUREMENTS 


5 


Proof. What we need to prove that for every H G 

E \{H,9x9l)Hsf<B\mls- 

AGZ2 

In other words, we need to prove that A||iJ||^g < ||yI(iJ)||^ < where A is the PhaseLift op¬ 

erator (I2.2I1 . Using the notation of the proof of Theorem 12.21 the formula (12.41) states that the fV-tuple 
S (C^)^ is obtained by performing inverse Fourier transforms of the columns of the matrix 
{N{H,g\gDqjg)x^^j^ = NA{H). Hence, their norms are related as follows: 

\\NA{H)f = ||(UX=ill' = ^ll(^’)f=iir- 


Using the same argument, we obtain 


(f"X=iir = ii(^^T=iir = 


-r,p\N 


”/V 


|K)^=ifand =N\\(nX=i\ 


The iV-tuples and are related through (12.51) . Therefore, if we define a = min;,gz 2 |(g,gA)|, 

P = max;^gz 2 |(g,gA)|, we have 

^ 2 

"'^nu\"= E E \Bpii){9,9p,)\ ^ 




> 




Finally, the matrix H is obtained by merely permuting the elements of the array Hence ||^f ||^5 = 

Combining everything, we obtain 


umx= 4 ii(^^)Eir=iiK),^=iir ^\ 

> 

which is exactly what we wanted to prove. 




= Np‘^\\{nX=i\^ =np‘^\\h 

= Na'^ \UnPXif =Na'^\\H 


\ls 

IL’ 


□ 


2.2.1. Non-vanishing vectors. We will now show that if we are interested in recovery of only non-vanishing 
vectors, weaker conditions on the generator can be assumed. 

Definition 2.2. A vector x G is called non-vanishing (or full), if all its entries are nonzero, i.e. 

x{n) 7 ^ 0, for all n = 0,..., fV — 1. 

By Cf we denote the class of all non-vanishing signals in C^. 

This situation is much easier to handle, because, intuitively, the non-presence of ’’holes” in the signals 
keeps the phases of the entries coupled. We will use the same technique as in Theorem l2.2l to prove that the 
injectivity condition can be weakened in this setting. Note that we are still assuming that all measurements 
are known. 

Theorem 2.4. Assume that 

(5; 9p,i) 7^ 0 for P = 0,1 and I G I^n- (2-6) 

Then the Gabor frame G = ig \)\^^2 allows phase retrieval for Cf. 

Proof. Assume that (12.61) is satisfied, and that x and y are full vectors which are measured equally by the 
Gabor frame. Then H := xx* — yy* is in the kernel of A, since for every A G Zfj- we have 

A{xx* - yy*)iX) = {gx, (xx* - yy*)gx) = |( 5 A,a:)|^ - \{gx,y)\^ = 0. 

The proof of Theorem 12.21 then implies that Tip = 0 for p = 0, 1, i.e. that TLo — TLi = 0. Remembering that 
TLp{i) = Hi^i-p, we arrive at 

0 = x{i)x{i) - y{i)y{i) = x{i)x(i - 1) - y{i)y{i - 1), i = 0,...,N -1. 
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The first equality simply says that |x(i)| = \y{i) \, i.e. that there exists numbers G T so that x{i) = eiy{i) 
for all i. Inserting this into the second equation yields 

0 = y{i)y{i - - 1). 

Since all entries of y are assumed to be nonzero, it follows that = e^-i, i.e. Ci = cq =: c G T for all i. 
Hence a: = cy for a c G T, and x and y are equal mod T. □ 

Remark 2.1. A similar result was proven in m- There, it was only assumed that {g, gp/) ^ 0 for p = 0 
and all £ G Zjv. However, in this case further constraints on the generators need to be made: g must be a 
window of length W >2, where N > 2W — 1 and N and W are coprime. Our result, on the other hand, 
works for more general generators and any N. 

2.3. Generators which allow phase retrieval. We will present two types of signals, one random and one 
deterministic, which satisfy condition (lOD . and thus can be used for phase retrieval of signals from all 
Gabor measurements. 


2.3.1. Complex random vectors as generators. We start by considering a probabilistic approach, a common 
strategy in signal recovery in general. 

Proposition 2.1. Let g be a vector in , randomly distributed according to the complex standard normal 
distribution. Then, the condition {g,gx) ^ 0 for all X G is satisfied with probability 1. 


Proof. Since there are only finitely many A’s, it suffices to prove that {g, n;,y) ^ 0 with probability 1 for one 
arbitrary A. Since H^ is a unitary operator, there exists an orthonormal basis of and G T with 

N 

Ha = '^c^q^q*. 

2 = 1 

If we expand g in this basis, i.e. g = hiqi, then the vector h G will also be distributed according to 
the complex standard normal distribution [T3]. We have Hap = Cihiqi, and hence 

N 

(5,nAy) = '^Ci \hi\‘^. 

2=1 


In order for g to not satisfy (ESI), the random variable h 
of K.^ defined by 



on must hence lie in the subspace 


V : 


ciVi = 0 

2 = 1 


Since this space has dimension N— 1, the set has Lebesgue measure zero. If we prove that t) has a distribution 
which has a density with respect to the Lebesgue measure on which is almost never zero, we are done. 
This is however not hard to see, since the variables \hi\^ = |ay + |&i|^, i = 1,...,N are independently 
distributed according to the y^'distribution, which has density p{x) = ^ exp{—x/2) on R+. □ 


2.3.2. Difference sets as generators. The second example are so-called difference sets, a construction coming 
from combinatorial design theory [S] . The set of all modulations of a characteristic function of difference set 
was shown to achieve the Welch bound in m- We will show that the set of all modulations and translations 
of a difference set has the property desired for phase retrieval. 

Definition 2.3. A subset K. = {ui,... ,uk} of is called an {N,K,v) difference set if the K{K — 1) 
differences 

(ufe — ui) mod N, k I 

take all possible nonzero values 1,2,..., — 1, with each value appearing exactly v times. 
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Example 2.1. Let N = 7. The subset /C = {1, 2,4} is then a (7, 3,1) difference set. We can check this by 
considering all possible differences modulo 7, 


- 

1 

2 

4 

1 

- 

6 

4 

2 

1 

- 

5 

4 

3 

2 

- 


and confirming that indeed every value from 1 to 6 appears exactly one time. 

Given a difference set /C with parameters {N, K, v) we denote by xk S {0,1}-^ its characteristic function: 


Xk{3) = 


ifj e/C 

if j^/C. 


We now prove that if such characteristic functions are used as generators, the corresponding Gabor frames 
will satisfy (lOD . and hence allow phase retrieval for arbitrary signals. 


Proposition 2.2. Let N he an integer with a prime factorization N = Let K. be a difference set 

with parameters (N, K, ly), such that 

u, iL < min{pi,...,pr}- (2.7) 

Then, for g = xk, 

7 ^^ for every (2.8) 

Proof. Let p, = {q,j), with both q,j 0. By just using the definition of g^ and K. we obtain 

{ 9 , 9 u) = ^ 9 (.n){MjTqg){n) = ^ g{n)g{n-q)u:^'^ = ^ (2.9) 

n^K. and 
n—q^K 

Now, taking into account the nature of a difference set, we can conclude that in the set 

{n : n € IC,n — q € K.} 

there will be always exactly v elements (because for g G Zjv there are exactly v ways to be written as a 
difference of elements in /C, and n — (n — q) are such differences). 

If u = 1, we are left with a single and then certainly the sum is different from zero. 

If u 7 ^ 1, we have a sum of v different iV-th roots of unity, and we will show that with the given assumptions 
on the difference set, (12.8|) holds. We use the following result from [TB] about the vanishing sums of roots 
of unity. The main theorem in this article states that for any N = pff .. .pff, the only possible amounts of 
iV-th roots of unity that can sum up to zero is given by Mipi +.. . + MrPr. Here the Mi are any non-negative 
integers (0 is included). Now it is clear that the condition v < min{pi ,.. .pr} will ensure that we will never 
have a vanishing sum. 

If now /i = (0, j) the sum will go over the full set /C, and since K < min{pi,.. .pr}, again this sum is non 
vanishing. 

Finally, in the last case p, = (^,0), we have a sum of v ones, and therefore we have proven (12.811 for all 
cases S □ 


Example 2.2. We now provide some examples of families of difference sets, which satisfy the condition 
from Proposition [221 

Family 1: Quadratic Difference Sets. Let q = p'' = Z (mod 4) be a power of a prime and 


N = q,K = 


g-1 

2 


u = 


g- 3 

4 


Then u = {ff : t G Zjv\{0}} is a {N, K, v) difference set. If r = 1, condition (12.711 is satisfied. 
Family 2: Quartic Difference Sets. Let p = 4a^ -I- I be a prime with a odd, and 


N =p,K = 


p-1 


p — 5 


4 


16 








PHASE RETRIEVAL FROM GABOR MEASUREMENTS 



1 

0 

1 



'1 

-1 






-1 


II 

1 



II 

0 
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Figure 1. The matrices Hi and H 2 used in the proof of Proposition [23] 


Then u = t G Z7v\{0}} is a {N, K, v) difference set and additionally K,i/ < N. 

Many other examples can be found in the paper or in the La Jolla Difference Set Repository at 
http;//www.ccrwest.org/ds.html, 

2.4. Generators which do not allow phase retrieval. We now consider two cases for which condition 
([231) is not satisfied, and show that this in fact implies that the Gabor frames do not allow phase retrieval 
in these cases. 

Proposition 2.3. Let g G he a generator such that one of the following two conditions is satisfied 

{9t9p,() = Oj for fixed p G ZAr\{0} and all £ G Z^r. (2-10) 

(ffj 9p,e) = 0: for p = 0 and all I G ZAr\{0}. (2-11) 

Then, the corresponding Gabor frame G = {g\)\^x‘^ does not allow phase retrieval for Cjv- 

Proof. Let us first assume that condition (I2.10p is satisfied. We consider the matrix Hi G defined by 

Hi = eoe*_p + C-pCp 

(eo is the ’first unit vector’ - remember that we are always considering indices from Zjv). This matrix has 
rank 2, and it lies in the kernel of the PhaseLift operator associated to the Gabor frame defined in (12.211 . To 
see this, note that using the notation of the proof of Theorem 12.21 we have 

fl iii = 0,p = p, 

TLpii) = Hip-p = % 1 if f = -p,p = -p, 

[0 else. 

In other words, Hp = 0 for all p ±p. Since {g,gp,e) = 0 J~^P{g,g-p-f), equation (12.101) also implies 
{g,g-p^e) = 0 for all £ G Z^v- These two facts prove that 

{9,9p,e) = 0 

for all £ and p. Using the technique of the proof of Theorem 12.21 backwards, it follows A{Hi) — 0. The 
matrix Hi that we have found has rank 2 and it is in the kernel of A. Therefore, by Theorem l2.ll the Gabor 
frame can not allow phase retrieval. 

Now we assume that (12.111) is satisfied. In this case we define a rank 2 matrix in by 

H 2 = eoCQ — eie*. 

For this matrix, Hp = 0 for p ^ 0. Also "££ 0 ( 0 ) = = 0. Because of these two facts and the assumption 

on g, we again have 

UpW (ff, 9p,e) = 0 for all (p, £) G Z^, 

and H 2 will by the same argument as before be in the kernel of A. Phase retrieval is again not possible. □ 
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Example 2.3. We now give two examples, for which the conditions of the previous proposition are satisfied. 

Short windows: The condition (12.101) is satisfied if the generator g is a “short window”. More precisely, 
if supp 5 C [Ki,K 2 \ for \Ki — K 2 \ < then g and M^Tpg will have disjoint supports for some p’s and 
hence have a vanishing scalar product. Using a window as a generator is a core idea in short-time Fourier 
analysis m- 


Alltop sequence: It can be easily shown that the much celebrated Alltop sequence [5], defined as 


^ )n=0 ’ 


has the property (12.111) . This generator is often and successfully used in sparse signal recovery from linear 
Gabor measurements I1I23]- 

However, both these families of signals can not be used for phase retrieval, when we are interested in 
recovery of all signals in C^. 


3. An injectivity condition for sparse signals 

3.1. Sparse phase retrieval via the PhaseLift operator. We will now consider signals which are sparse 
in a dictionary. A dictionary D is & set of d vectors in C^, and it is identified with the matrix formed when 
writing the d vectors as its columns. The class of signals which are fc-sparse in the dictionary H, or simply 
fcH-sparse signals, is 

Ck,D = {x S C^l 3 z € C^, ll^llp < k, s.t. X = Dz] , 

where ||2:|1 q denotes the number of non-zero coefficients in z. If H = J, we will omit the dictionary and 
simply speak of fc-sparse vectors Ck- We will also, to increase readability, speak of fcZ?-phase retrieval instead 
of phase retrieval for Ck,D- 

Since sparse vectors in some sense are k- and not A^-dimensional, one would hope that the number of 
measurements required to recover them is smaller (for us, A should contain less elements.) This, and other, 
questions were considered and answered for general measurement vectors in [22) and [26]. A counterpart 
of Theorem [Q in the sparse setting has up to know not been stated and proved. We will prove now an 
injectivity condition for sparse signals, and then use it for the case of Gabor measurements as in the previous 
section. 

For a given dictionary D = {di)‘^^^ and a set of indices /C C {1,... d}, we denote W/c = span{di}ig;c- With 
the help of this notion, we can characterize sets or measurement vectors which allow fcH-phase retrieval. Let 
A be the PhaseLift operator defined in (12.21) . 

Theorem 3.1. Given the notations from above, the following two statements hold. 

(1) If for every K, with |/C| = 2k, the kernel of A does not contain rank 1 or 2 Hermitian matrices whose 

range is in Wic, then the vectors allow kD-phase retrieval. 

(2) If (/i )™2 allowing kD-phase retrieval, then for every K, with \K\ = k, the kernel of A does not 
contain rank 1 or 2 Hermitian matrices with range in W)c. 

Proof. Let us start by proving (1) by contraposition. Assume that (/i)™ is not allowing fcH-phase retrieval. 
Then there exists x y mod T, both kD-sparse, for which 

{xx*JJ*)hs = \{fi,x)f = \{fi,y)\^ = {yy*Jift)Hs^ 

i.e. XX* — yy* is a Hermitian matrix in the kernel of A. If the sparse representations of x and y are given by 
X = Dzx and y = Dzy, we see that Ya\i(xx* — yy*) C (Wsuppzj,usupp 2 „) and further it has rank less than or 
equal to 2. If we knew that the rank is at least one, |supp Zx U supp Zy \ < 2k would imply the claim. 

To see that the condition x ^ y mod T in fact implies this, assume, towards a contradiction, that this is 
not the case, i.e that xx* — yy* = 0. Since x ^ y mod T, both vectors are non-zero. Hence there exists a 
vector V G such that {x, v) 0. Multiplying 0 = xx* — yy* with this v and rearranging terms, we arrive 
at 

{y,v) 

X = 7 - ry, 

{x,v) 

i.e x = Ay for a A G C. Again plugging this into 0 = xx* — yy* yields |A| = 1. This is a contradiction. 
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Let US now turn to (2). Suppose that there exists a /C such that the kernel A contains a Hermitian matrix 
H with rank 1 or 2 with range in W/c- By the spectral theorem, there exists an orthonormal basis {ipj) of 
consisting of eigenvectors of H, corresponding to real eigenvalues (A^). It is clear that H may be written 
as 

Because of the bounded rank and the fact that H 0, either one or two of the eigenvalues are non-zero. 
It is clear that the eigenvectors corresponding to those eigenvalues are vectors in Wic , since they form a basis 
of the range of H. Thus, they are fciA-sparse. 

Let us first consider the case where only one eigenvalue is different from zero. If we write x = 
then we have H = ±xx* and hence 

0 = (xx*Jif*} = \{fi,x)f = 0. 

This means that the two fcZI-sparse vectors x and 0 have the same phaseless measurements, although x ^ 0 
mod T. 

The other case is dealt with similarly: here we write x = y = and conclude that 

H = ±xx* ± yy*, where the signs depend on the signs of the eigenvalues. If the signs are equal, we see that 
|(/i,x)|^ + |(/iiy)l^ = 0 and we have again found fcD-sparse vectors which are measured 0. If the signs are 
not equal, we see that \{fi,x)\ —\{fi,y)\ = 0, and hence x and y are measured equally. They are fciA-sparse 
and cannot be equal mod T, since they are orthogonal. □ 

3.2. Signals sparse in the standard basis. Let us start by considering vectors which are sparse in the 
standard sense, i.e. D = I. We will prove a condition under which a subset of our Gabor frame {g\ , A G } 
with ^ elements allows fc-sparse phase retrieval, when N is prime. For general N, it would still be possible 
to go below the full set of measurements, N'^. We will need a special form of the discrete uncertainty principle, 
which involves the sum of the “spread” of the signal and its Fourier transform. Let us start with a general 
observation. 

Lemma 3.1. Assume that for all non-zero vectors f G 

||/|lo + ll/llo>A^-0iv (3.1) 

holds for some number 0^. Then, if f is k-sparse (l|/||o = k), and / has not less than 0N + k-\-l zero-entries, 
then f necessarily has to vanish. 

This statement follows immediately by contradiction. The question is whether en) is a reasonable 
assumption. In |25j it is proved that when N is prime, ii]) holds with 9n equal to —1. For general N, by 
the standard multiplicative uncertainty principle and the geometric mean-arithmetic mean inequality, one 
can derive (13.11) with 9^ = N — 2^/N. A more involved inequality for general N was obtained in m and 
will be discussed later on. 

Before we proceed with a condition on the generator g for sparse phase retrieval, we will first prove a 
more general statement about recovery of sparse matrices from linear measurements, which is interesting on 
its own. We will be interested in the following class of signals, 

S^K = {iL G : 3/C C [l,...iV],|/C| = K : H,j = 0 if {i,j) ^ /C x /C} . 

Theorem 3.2. Let N be such that the uncertainty principle (ED holds, and let A = {p, 1) G Z^. Let g have 
the following property: for each £, the sequence Cp = ((ff, ffp,-)) formed by letting I run obeys 

9N + K + \<\\cp\\^<k (3.2) 

for some K and k. Then, for any subsets BCZjv with 

|A| > 6»A,-kfc-kl, \B\>9n + K'^-K + 2, 

the following holds. If a matrix H G satisfies {{gxg’^, H)^g)\^AxB = 0, then H = 0. 
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Proof. Let H G satisfy {gxg^, H)^g = 0 for A G x B and let fC be such that Hij = 0 if {i,j) ^ 1C x 1C. 
We will prove that H then must be 0. Recall the notation from the proof in Theorem 12.11 ’Hp(i) = 

Since Hi^i-p is zero, if (i, i — p) is not in /C x /C, we can conclude that 

np{i) = H,^,_p = Q ifi ^/Cn(/C+p). 

This proves the following properties: 

(1) The vectors Up are RT-sparse. 

(2) Up is zero for all but at most — K + 1 different values for p. To see this, notice first that Up = 0 

if p ^ K. — 1C. This is because if 'Hp(i) ^ 0, then i ^ K. and there additionally exists a j G /C 
with i = j + p. It follows p = i — j G K. — 1C. And we know that the set K. — K, has at most 
|/C| (|/C| — 1) + 1 = — AT + 1 elements. 

Now using the same argument as in the proof of Theorem 12.21 we arrive at 

0 = N{g^glH)=Y,^~^^^''^npi£){g,gp,i) for all A = (g,j) G A x B. (3.3) 

p,£ 

Fixing J, the sum in (liTATl) is the value at j of the discrete Fourier transform of the vector defined as 

V^{p)=Y,Jmp{e){g,gp,f). (3.4) 

i 

Because of these vectors are all (AT^ — AT + l)-sparse. Further, (13.31) proves that their Fourier transforms 
vanish at all j G B, i.e. at 6(jv + (AT^ — AT + 2) points. The discrete uncertainty principle (13.11) implies that 
must equal zero. 

Considering (El, the fact that V^{p) = 0 proves that the inverse Fourier transform of the vector, which 
we denote by 

wP{e) = Upit) {g,gp,e) 

vanishes at the values q € A, i.e. at + A: + 1 values. Because of our assumption on p, is however 
^-sparse. We can therefore again conclude that 

{g,gp,e) = 0 for all (p,£) G 

Hence, if {g,gp,i) ^ 0, 'Hp{t) must be 0. Due to our assumption on g, this happens for at least + 2fc + 1 

fs for every p. Because of[Tl this is sufficient to prove that Up = 0 for all p, and H therefore must be 0. □ 

We now use the theorem we have just proved, to provide a condition, when a Gabor frame can do fc-sparse 
phase retrieval. 

Theorem 3.3. Let N be such that the uncertainty principle (EH) holds, and let A = (p, 1) G Z^. Let g G 
be a generator which satisfies the following condition: for each I, the sequence Cp = ((p,Pp,.)) formed by 
letting £ run obeys 

9N + K + l<\\cp\\Q<k (3.5) 

for some K = 2k and some k. Then, for any subsets A G 'Z]\j, B Q Zat with 

|A| > 0Af + fc + 1, \B\ > 0Ar + (2A:)^ — 2k + 2, 

the set 

{gx, XcAxB} (3.6) 


allows k-sparse phase retrieval. 
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Proof. We will use part (1) of Theorem 13.II for D = I, to show that fc-sparse phase retrieval is possible for 
the system (EH). Let H be an Hermitian operator with values in = {x G , supp(a;) C jC} for some /C 
with |/C| = 2k for which A{H) = 0 (where A is the PhaseLift operator associated with (I3.6p l. We will prove 
that H must be zero, from which the claim follows. Since the range of H is contained in C^, we know that 
Hij = 0 if i ^ /C. Since Hij = Hj^i, we also have Hij = 0 if j ^ /C. We can conclude that H G ^ 2 K, and 
by Theorem 13.21 it immediately follows that H is zero, thus the theorem is proved. □ 

Remark 3.1. We would like to state the following remarks related to Theorems 13.31 and 13.21 

(1) If 0Ar + 4fc^ — 2fc + 2 > N, then the same theorem holds for B = Z^r and any A with |a1| > 9n + k + l. 

We can therefore also in this case reduce the number of measurements from to {9n + k+l)N. Also 

note that since k must not be smaller than 2k, the theorem does not yield any enhanced results for 
non-sparse vectors (we need the sequences Cp to be A:-sparse to reduce the number of measurements, 
but to have at least 2k nonzero elements to ensure injectivity for fc-sparse signals). 

(2) We note, that when N is prime, the conditions of Theorem 13.31 become much simpler. Namely, 

2fc < llcpilo < fc, 

and the sets A and B should fulfill 

|2l|>fc, \B\ > {2kf - 2k + 1. 

Thus, for example, if we can find a Gabor system for which the inequality (1531) is fulfilled as an 
equality, we will be able to do fc-sparse phase retrieval with order of only 0(fc min(iV, fc^)) measure¬ 
ments. 

(3) When N is not prime, we have 9^ = N — 2\fN, and the number of needed measurements is not as 
good as in the prime case, since we obtain 

|A| • \B\ > {N - 2^/N + k + 1){N - 2^/N + {2kf - 2k + 2), 

but some improvement over could still be obtained in some cases. Furthermore, an extension of 
|25j from N prime to general N was published in |19) . in the form of the following property: 

Let di < 9,2 be two consecutive divisors of TV. If di < fc = ||/||g < 9,2, then 

ll/llo > +d 2 - fc). 

did2 

Our function 9 will in this case explicitly depend on fc and be equal to N — k + j^idi + d 2 — k). 
The smaller this value is, the less measurements will be needed for fc-sparse injectivity. 

(4) Theorem 13.21 is interesting from a different perspective, since can be viewed as a set of K^- 
sparse vectors in C whose sparsity has a special structure. Thus, we have provided a deterministic 
construction which can theoretically recover those vectors from O(K^) linear measurements. This 
is interesting since we know from conventional compressed sensing results |12] . that deterministic 
constructions for stable recovery of iV^-sparse vectors require 0{K^) linear measurements, whereas 
random constructions only need 0{K^) measurements. Finding deterministic constructions which 
can accept sparsity levels on the order higher than the square root of the number of measurements 
is known as breaking the “square-root bottleneck” [50] . 

Although in our case, 0{m) measurements are needed for sparsity level rnfA^ one has to bear in 
mind that the sparsity of the vectors is structured, and that our result is only about the injectivity of 
the measurements. In particular, we do not prove any recovery guarantees for a specific algorithm. 
Hence, we have not broken the square-root bottleneck, but the theorem can be seen as a step towards 
providing new results in this direction. 

3.3. Functions window in the Fourier domain as generators. As in the previous section, we now 
provide an example of a generator g which fulfills our condition. 
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Proposition 3.1. Let N be prime and 2k + 1 < N. Further, let v G be a window of length fc + 1, 
V = X[o,k\! where XA denotes the eharacteristic function on the set A C Zjv- Moreover, let g be defined by 
g = V, Then, g satisfies (ixa with k = 2k + 1 and therefore, (2k + 1) inin(4A;^ — 2k + 1, N) measurements 
from the Gabor frame with g as generator will do k-sparse phase retrieval. 


Proof. The Plancherel formula implies that 

Cp(^) = {g,TiMpg) for all {p,t) G 


Therefore, 

(5,r^Mpg) = ^ ^ ^jP^rn-e)^ 

m^[0,k] and 
m—t^ [0,fc] 


(3.7) 


because of the way we defined v. Note that since 2k < N, this sum is empty for \i\ > k. Therefore, the 
sequences Cp are 2k + 1-sparse for every p, i.e. A:-sparse. 

It remains to prove that for \£\ < k, the expression above is not zero, and hence ||cp||p = 2k + l. It suffices 
to consider .^ > 0, since the other case can be obtained from this one by the substitution i -A —1. Using the 
formula for geometric sums, we obtain 


£<m<k 


k-e 

n—0 


1 _ ijjp(k-e+i) 

1-UjP 

k — £ + 1, 


iipy£0, 


if p = 0. 


The only way this could be zero when £ < k is that p 0 and 1 — = 0. This would however mean 

that TV is a divisor of p{k — £ + 1) ^0. Since N is prime and both p and (fc — £ + 1) are smaller than N, this 
cannot be the case. Therefore, from Theorem liOl we conclude that any subsets ACZjv, SCZ^v with 

\A\>2k + l, \B\ > {2kf - 2k + 2 


will do fc-sparse phase retrieval. 


□ 


Remark 3.2. The choice of p as a characteristic function of [0, k] is not necessary - any generically chosen 
function with support on [0, k] will also lead to a Gabor system with the same properties. To see this, note 
that if |£| < k, the expression (13.71) is a non-trivial polynomial in the variables re(w),im(w). Since we have 
proved that there is a particular choice of v so that all polynomials do not vanish on v, (13.51) will be satisfied 
for generic v. 

Remark 3.3. The matrices provided to prove that the frames considered in Section [231 do not allow phase 
retrieval were all matrices with range in for a /C with |/C| = 2. Hence, the considerations made there in 
fact proved that the frames are not allowing phase retrieval for Ck for any k > 2 (although they might still 
allow phase retrieval for some other class of signals). 

3.4. Signals sparse in Fourier domain. After spending some time discussing the standard sparsity case, 
it is worth noting that similar results hold for signals which are sparse in the Fourier basis (dictionary) F. 
Recall the famous commutation relation of F with translations and modulations: 

n(p,f)F = MeTpF = FTeM.p = w^PFn(_p,^) for all (p, £) S (3.8) 

This formula allows us to translate the results provided in the previous section to this new setting. 

Theorem 3.4. Let N be such that the uncertainty principle (EH) holds and F denote the Fourier basis. Let 
g have the following property: for each £, the sequence eg, = {{g,g.,g)) formed by letting p run obeys 

9]\[ + 2fc + 1 < ||cf llg < k (3-9) 

for some k and k. Then, for any subsets A, B Q Z^v with 

|A| > 9i\[ + (2fc)^ — (2£:) + 2, |R| > 9]\[ + fc + 1, 

the set {g\, A S A x B} allows Fk-sparse phase retrieval. 
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Proof. We would like to apply Theorem 13.11 with D = F. Using the notation of that theorem, let H be 
an arbitrary Hermitian matrix with range contained in W/c for some K, with |/C| = 2k. We may write 
H = FH^F* for some other Hermitian which then has a range which is contained in C^. Let us now 
proceed as in the proof of Theorem 12.21 and calculate 

= = co-^pn^M- 

We used the commutation relation (EHl), and the fact that F is unitary. We arrive at 

N {g\, Hgx) = {gx, Ii(p,i)g\) ■ 

p,i 

This formula is very similar to (13.31) . essentially, the only difference is that the roles of p and i have inter¬ 
changed. Further, the vectors have the same properties as the vectors PL in the proof of Theorem 13.41 
(since has the same properties a.s H). These two facts makes it clear that we can use the exact same 
technique as in that proof to prove this theorem. We leave the details to the reader. □ 

It is not hard to construct a concrete example of a generator g which fulfills the condition (13.91) . We only 
have to note that the roles of translations and modulations have been interchanged. Hence, we should no 
longer use a g which has short support in Fourier domain, but instead one with short support in spatial 
domain. With this insight, we may use the exact same steps as in the proof of Proposition [3T] to deduce the 
following. 

Proposition 3.2. Let N be prime and 2k + 1 < N. Then {2k + 1) min(4A:^ — 2A: -I- 1, N) measurements from 
a Gabor frame generated by generic windows g of length k + 1 will do Fk-sparse phase retrieval. 

4. An ALGORITHM FOR PHASE RETRIEVAL USING GABOR MEASUREMENTS 

The idea of the proof of Theorem l3.3l can be used to design an algorithm to reconstruct signals from their 
Gabor intensity measurements. We start by recovering F[, as in the proof, and then we compute the closest 
rank 1 operator xx* by spectral decomposition of FI. A detailed description is given in Algorithm 1. 


Algorithm 1: Simple Gabor Phase Retrieval (SGPR) 

Data: A generator g G C^, sets A,B C Z^r, the measurements b{q,j) = N\{x,gqj)\'^ , {q,j) G Ax B. 
Result: An estimate xq G of x. 

1 for q = 0 ... N — 1 do 

|_ Solve U«(j) = b{q,j), j G B for U®. 

2 for p = 0 ... A — 1 do 

|_ Solve N ■ wP{q) = U‘^(p), qG Aior w^. 

3 for p = 0 ... A — 1 do 

for f = 0 ... A — 1 do 

if {g, 9p,i) ^ 0 then 

Set 'Hp(f) = wP{l)l {g, g^^i). 

Add t to the set Ap. 

Solve 'Hp(f) = wP{i)f {g,gp,i), i G Ap for Hp. 

Reconstruct H from H{i,i — p) = T~Lp{i) 

Calculate the eigenpair (A, v) of FL corresponding to the largest eigenvalue. 

Set Xq = y/Xx. 


In steps (I), (2) and (3) one has to invert a Fourier transform. If all values of the transformed vector are 
known, one can simply use the standard fast inverse Fourier transform to compute this, and the signal will be 
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perfectly recovered. If one on the other hand does not know all the values (not all Gabor measurements are 
given), some other method has to be used, where sparsity can be employed. We have chosen Basis Pursuit 
[5] . This is a standard approach in compressed sensing when looking for a sparse solution x of the equation 
Ax = h. The algorithm consists of solving the following optimization problem: 

min ||x|| s.t. Ax = b. 

We solve this problem with CVX, a package for specifying and solving convex programs M- All experiments 
were conducted on a Intel Core i7-3517U Processor running Windows 8. 



o'-'-^-'-'-^-'-' 

1 2 3 4 5 6 7 8 

Sparsity level, k 


(A) Partial modulations only 


(B) Partial modulations and translations 


Figure 2. Success rate 


We now present the results of the numerical experiments for testing Algorithm 1. In Figure [21 we plot 
the success rate of recovery of sparse signals via Algorithm 1. We have fixed the length of the signal N = 67 
(a prime which gives 3 mod 4, as needed for difference sets of Family 1). We want to recover two types of 
signals: fc-sparse signals, where k non-zero random values are distributed on a random support, and A:-sparse 
block signals, where one random block of k subsequent entries is assigned k random values. 

In Figure [2^ we also have chosen two different generators for the Gabor system: a complex random signal, 
and a characteristic function of a difference set, described in Section |2l Here, we use 0.5fV^ measurements, 
namely all N translations, and random 0.57V from the modulations. With this setup, we use ii minimization 
only in the Step (1), and in (2) and (3) we use the fast Fourier transform. For a fixed sparsity from 1 to 15, 
we repeat the experiment T = 200 times, and count a trial as successful, if the normalized squared error is 
smaller than 10“^. 

In Figure [2Bl we do the same experiment, but we take partial measurements in both directions: translation 
and modulation. Namely, for N = 67, we take 0.52A^ translations, and 0.7iV modulations at random. The 
generator here, as described in Proposition 13.11 is a short Fourier window, with length 8. Now, we need to 
use Basis Pursuit in all steps (1), (2) and (3), which in turn leads to a lower recovery rate. We made T = 100 
trials for every sparsity level. 

In Figure [351 we test the speed of our algorithm in comparison to the PhaseLift algorithm [B], implemented 
using the CVX package. We also use Gabor measurements for it, but only 21og(7V)7V, taken at random. 
We plot the average execution time over T = 50 trials, and see that as the dimension grows, our method 
becomes faster, although the number of measurements is much larger. Also, if we are using the full set of 
measurements, the time needed is incomparably smaller to both of the other methods - since then there is 
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(A) SGPR VS. PhaseLift 


(B) SGPR vs. GESPAR 


Figure 3. Recovery time 


no minimization problem included. In this case, also, we will always recover the signal with probability 1, 
independently of the sparsity level. 

In Figure l3Bl we compare the execution time of Algorithm 1 from all measurements to the GESPAR 
algorithm m. a greedy algorithm for recovery of sparse signals from Fourier intensity measurements (in 
our experiment we use 2N measurements). This algorithm is very fast for high dimensions, but since it is 
iterative, it becomes slower as the sparsity increases for a fixed dimension of the signal. We illustrate this 
behavior in Figure [SB] where for every dimension, we measure the average time of recovery of signals with 
sparsity fc = 5 and k = 10. We see that the GESPAR algorithm is faster, when we want to recover a signal 
which has only 5 nonzeros, but if this number is larger, our algorithm becomes faster than the GESPAR, 
since it does not strongly depend on the sparsity level. 

We would like to mention that our algorithm for all iV^ measurements is also stable to additive noise in 
the measurements. This follows from Theorem l2.3l and can be intuitively explained by the fact that the only 
troublesome part is the division in Step (3). If the generator g is such, that the values {g,gp,i) are bounded 
away from zero, one can guarantee robustness to noise. For the recovery from less than measurements, 
we leave the detailed investigation on this matter for future work. 
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